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SUMMARY 


The  purpose  of  this  work  was  to  develop  and  evaluate  an 
interactive  long-period  signal  extraction  filter  suite.  Five 
filters  were  included  in  the  suite  and  were  applied  to  test 
data  in  a  cascaded  mode,  (i.e.,  the  output  of  one  filter 
served  as  the  input  to  the  next  filter).  The  filters  tested 
were  a  bandpass  filter,  a  Wiener  filter,  a  chirp  matched 
filter,  a  polarization  filter,  and  a  Gaussian  narrowband 
filter.  It  was  found  that  for  the  filters  examined,  the 
polarization  filter  dominates  the  cascaded  filtering  approach 
to  signal  extraction;  specifically,  processing  gains  as  high 
as  14  to  IS  dB  were  realized  from  the  application  of  cascaded 
filters  which  included  this  filter.  Further,  tests  performed 
on  a  Eurasian  data  set  using  one  specific  cascaded  filter 
sequence  yielded  an  improvement  in  surface  wave  detection 
capability  equivalent  to  an  improvement  in  bodywave  detec¬ 
tability  of  approximately  0.6  to  0.8  mj^  units.  Finally, 
cascaded  filters  were  found  to  provide  more  detected  signals 
from  which  surface  wave  magnitudes  can  be  obtained;  as  such, 
the  use  of  these  filters  can  enhance  an  analyst's  capability 
to  discriminate  between  earthquakes  and  explosions. 

Neither  the  Advanced  Research  Projects  Agency  nor  the 
Air  Force  Technical  Applications  Center  will  be  responsible 
for  information  contained  herein  which  has  been  supplied  by 
other  organizations  or  contractors,  and  this  document  is 
subject  to  later  revision  as  may  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should 
not  be  interpreted  as  necessarily  representing  the  official 
policies,  either  expressed  or  implied,  of  the  Advanced  Re¬ 
search  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  US  Government. 
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SECTION  I 
INTRODUCTION 


A.  THE  TASK 

This  report  presents  the  results  o£  a  study  performed 
on  cascaded,  long-period  signal  extraction  techniques.  The 
purpose  of  the  study  was  fourfold.  First,  it  was  desired  to 
integrate  a  suite  of  long-period  filters  into  an  interactive 
seismic  processing  system.  Second,  the  procedures  for  using 
the  filters  in  a  cascaded  mode  were  to  be  standardized  and 
optimized  for  automatic  mode  operation.  Third,  the  filter 
program  was  to  be  tested  for  its  effect  on  signal  detection 
and  event  discrimination.  Fourth,  the  feasibility  and  oper¬ 
ability  of  using  the  long-period  signal  extraction  program 
on  a  large,  world-wide  data  base  were  to  be  determined. 

Early  in  the  course  of  the  study,  a  problem  arose  which 
forced  a  change  in  the  manner  in  which  the  work  was  performed. 
Specifically,  the  problem  was  that  computation  of  a  fast 
Fourier  transform  (FFT)  on  the  PDF- 15/50  (an  interactive 
processing  system)  is  significantly  slower  than  it  is  on  the 
IBM  360/44  (batch  processing  system).  Since  the  polarization 
filter  used  here  requires  many  FFTs  in  order  to  compute  and 
apply  the  appropriate  filter  weights,  the  relatively  slow 
computational  speed  of  the  interactive  system  rendered 
impractical  the  processing  of  large  data  bases. 


ENSCO,  INC. 


I-l 


To  circumvent  this  problem,  the  interactive  filter  pro¬ 
gram  was  primarily  used  to  estimate  the  optimum  values  of 
these  filter  parameters  which  had  not  been  determined  in 
earlier  signal  extraction  studies  (Strauss,  1978).  The  in¬ 
teractive  system  was  also  used  to  test  the  filter  algorithms 
as  each  was  designed.  However,  to  perform  tests  of  the  cas¬ 
caded  filters  on  large  data  bases  a  similar  filter  program 
suite  was  developed  for  use  on  the  batch  processing  system 
(hereafter  referred  to  as  the  automatic  mode  of  processing) 
using  the  optimum  values  of  the  controlling  parameters  which 
were  determined  by  means  of  the  interactive  system. 

It  should  be  noted  that  the  problem  described  above  (the 
FFT  computational  speed  on  the  PDP-15/50)  is  not  common  to 
all  interactive  systems.  A  transfer  of  the  programs  develop¬ 
ed  here  to  a  system  such  as  a  PDP-11/70  with  an  array- 
processor  subsystem  would  provide  for  efficient,  interactive 
processing  of  large  data  bases.  The  transfer  of  the  inter¬ 
active  signal  extraction  program  to  another  system  remains 
as  a  subject  for  future  study. 


B.  SUM^fARY  OF  PREVIOUS  WORK 

A  number  of  studies  have  been  made  on  the  use  of  specific 
filtering  techniques  for  long-period  signal  extraction.  Rele¬ 
vant  results  obtained  in  some  of  these  studies  are  summarized 
in  Table  I-l.  In  each  case,  the  signal  extraction  technique 
was  applied  to  'raw'  (i.e.,  previously  unfiltered)  data.  Note 
that  of  the  filters  examined,  the  greatest  improvement  in  the 
signal -to-noise  ratio  (SNR)  was  obtained  using  the  polariza¬ 
tion  filter  (7.5  dB  to  10  dB  for  Rayleigh  waves). 
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TABLE  I-l 


SUMMARY  OF  PREVIOUS  SIGNAL  EXTRACTION  RESULTS 
USING  VARIOUS  FILTERING  TECHNIQUES 


Extraction  Technique 

Process¬ 
ing  Gain 
rdB'l 

Source  Of 
Test  Data 

Reference 

Chirp  Matched  Filter 

3 

Eurasian 
earthquakes 
recorded  at 
a  single 
site  of  the 
Alaskan 

Long  Period 
Array 

Strauss,  1973 

Reference -Wave form 
Matched  Filter 

3 

Strauss,  1973 

Polarization  Filter: 
Rayleigh  Waves 
Love  Waves 

7.5 

6.5 

Strauss,  1976 

Wiener  Filter 

2-6 

Kurile 

Islands  - 
Kamchatka 
earthquakes 
recorded  at 
the  Guam 
Seismic 
Research 
Observatory 

Lane,  1976 

Bandpass  Filter  - 
Polarization  Filter 
(Rayleigh  Waves) 

10 

Lane,  19  77 

Wiener  Filter  - 
Polarization  Filter 
(Rayleigh  Waves) 

10 

Lane,  1977 

Wiener  Filter  - 
Polarization  Filter 
(Rayleigh  Waves)- 
Re ference -Wave form 
Matched  Filter 

8 

Lane,  1977 
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The  first  results  obtained  using  cascaded  filters  of  the 
types  discussed  here  (Lane,  1977)  are  also  presented  in  Table 
I-l.  Note  that  when  a  reference -wave form  matched  filter  was 
cascaded  with  other  extraction  techniques  it  lowered  the  net 
gain  of  the  cascaded  filters  by  2  dB.  This  result  may  be 
due  to  the  difficulty  encountered  in  selecting  a  suitable 
reference  event  from  a  complex  source  region  such  as  the 
Kurile  Islands.  Based  on  this  result,  chirp  matched  filters 
were  substituted  for  reference -wave form  matched  filters  in 
the  present  study. 


C.  SPECIFIC  GOALS 

The  specific  goals  of  this  study  are  to: 

•  Integrate  available  long-period,  signal-extraction 
filtering  programs  into  an  interactive  seismic 
processing  system; 

•  Standardize  and  optimize  the  procedures  for  using 
the  interactive  seismic  processing  system; 

•  Create  an  automatic  mode  processing  version  of  the 
interactive  program; 

•  Estimate  SNR  improvement  produced  by  the  applica¬ 
tion  of  cascaded  filtering  techniques; 

•  Estimate  the  effects  of  cascaded  filtering  tech¬ 
niques  on  detection  capability; 

•  Estimate  the  effects  of  cascaded  filtering  tech¬ 
niques  on  earthquake/explosion  discrimination 
capability ; 
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•  Determine  the  feasiblity  and  operability  of  using 
the  new  long-period  filtering  program  on  a  large, 
world-wide  data  base. 

The  work  performed  to  fulfill  each  of  these  goals  is 
presented  as  follows.  Section  II  presents  a  description  of 
each  of  the  filtering  techniques  used  and  of  the  method 
used  to  combine  these  techniques  into  interactive  and  auto¬ 
matic  mode  programs.  Also  included  in  this  section  is  a 
description  of  the  controlling  parameters  peculiar  to  each 
filtering  technique. 

Section  III  presents  information  on  the  SNR  improvement 
obtained  from  various  sequences  of  cascaded  filtering.  A 
discussion  of  the  optimum  cascaded  filter  sequence  is  also 
presented. 

Sections  IV  and  V  present  the  results  of  an  experiment 
performed  on  a  large  data  base  using  one  specific  sequence 
of  cascaded  filters.  The  results  of  the  experiment  illus¬ 
trate  the  effects  on  signal  detection  and  discrimination 
produced  by  cascade  processing.  Specifically,  Section  IV 
discusses  the  improvement  obtained  in  the  detectability  of 
seismic  signals  recorded  by  single  stations  and  by  a  small 
network  of  stations.  Section  V  discusses  the  effect  of  the 
filter  sequence  examined  on  earthquake/explosion  discrimina¬ 
tion  by  means  of  the  M^-m^  discriminant. 

Finally,  Section  VI  summarizes  the  results  of  this 
study  and,  based  on  the  results  obtained,  discusses  the 
feasibility  of  using  cascaded  filters  on  a  large,  world¬ 
wide  data  base. 
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D.  TEST  DATA  BASE 


Test  events  used  in  studies  such  as  this  should  be 
located  in  the  same  region  in  order  to  minimize  the  effects 
of  source- receiver  path  differences  on  event  detection  and 
discrimination.  Accordingly,  approximately  200  earthquakes 
which  occurred  between  August  1978,  and  April  1979,  and  which 
had  epicenters  in  the  Kashmir-Uzbek-Tadzhik  area  (as  listed 
in  the  Norwegian  Seismic  Array  event  bulletin)  were  selected 
for  the  earthquake  test  data  base.  Since  no  region  yielded 
an  appreciable  number  of  presumed  explosions  during  the  time 
frame  noted  above,  all  Eurasian  presumed  explosions  were 
selected  for  use  in  the  discrimination  test.  Data  for  all 
events  were  edited  from  the  digital  magnetic  tapes  recorded 
at  the  Ankara,  Turkey  (ANTO) ;  Shillong,  India  (SHIO) ;  Taipei, 
Taiwan  (TATO) ;  Matsushiro,  Japan  (MAJO) ;  Guam,  Marianas 
Islands  (GUMO);  Chiang  Mai,  Thailand  (CHTO) ;  Kongsberg,  Nor¬ 
way  (KONO) ;  and  Grafenburg,  Germany  (GRFO)  Seismic  Research 
Observatories.  The  source  regions  and  recording  stations 
are  shown  in  Figure  I-l,  where  a  solid  circle  (•)  represents 
a  station,  a  star  (•)  represents  a  presumed  nuclear  explosion 
(PNE)  source  location,  and  the  solid  square  (■)  represents 
the  earthquake  source  region.  The  PNE  locations  were  taken 
from  bulletins  of  the  National  Earthquake  Information  Service 
(NEIS)  . 
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SECTION  II 

THE  LONG -PERIOD  SIGNAL  EXTRACTION  PROGRAM 

/ 


A.  GENERALIZED  STRUCTURE  OF  THE  PROGRAM 

The  general  structure  of  the  long-period  signal  extrac¬ 
tion  program  is  shown  in  Figure  II- 1.  Both  the  interactive 
and  automatic  mode  versions  can  be  considered  to  have  four 
processing  stages:  (1)  prepare  and  display  data;  (2)  select 
filter;  (3)  apply  filter;  and  (4)  display  filtered  data  and 
make  magnitude  measurements.  Both  versions  were  designed 
for  ease  of  insertion  of  additional  signal  extraction  algo¬ 
rithms,  For  the  interactive  program,  this  is  accomplished 
through  user  functions  (Shaub  and  Black,  1977)  while  for  the 
automatic  mode  program,  this  is  accomplished  through  the  use 
of  subroutines. 

Because  of  the  FFT  time  constraint  discussed  in  Section 
I,  the  interactive  system  was  used  primarily  to  determine  the 
optimum  manner  in  which  to  apply  the  automatic  mode  version 
of  the  program.  Thus,  each  processing  stage  of  the  interac¬ 
tive  program  shown  in  Figure  II- 1  requires  analyst  interven¬ 
tion  in  the  form  of  time  gate  and  filter  parameter  selection. 
The  automatic  mode  version  of  the  program,  on  the  other  hand, 
requires  no  such  intervention.  Processing  time  gates  are 
selected  on  the  basis  of  expected  waveform  arrival  times 
while  filter  parameters  are  preset. 
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Stage  3 


Stage  4 


FIGURE  II-l 

PROGRAM  STRUCTURE  FOR  CASCADED  FILTERING 


B.  SIGNAL  EXTRACTION  TECHNIQUES 


This  section  presents  brief  descriptions  of  the  five 
filter  techniques  used  in  this  study  and  the  manner  in  which 
they  are  combined  to  form  the  long-period  signal  extraction 
program.  The  five  filter  techniques  used  are: 

•  Bandpass  filter  with  cosine-squared  taper 

•  Wiener  filter 

•  Chirp  matched  filter 

•  Polarization  filter 

•  Gaussian  narrowband  filter. 

1.  Bandpass  Filter  with  Cos ine- Squared  Taper 

The  bandpass  filter  operates  in  the  frequency  domain, 
suppressing  all  power  at  frequencies  below  the  low  bandpass 
frequency  (BPL)  and  above  the  high  bandpass  frequency  (BPH) 
(see  Figure  II-2).  The  algorithm  applies  a  cosine-squared 
taper  to  the  power  at  those  frequencies  between  BPL  and  BPL- 
plus-the-taper-length,  and  between  BPH-minum-the-taper- length 
and  BPH.  (The  taper,  which  is  used  to  minimize  filter  ring¬ 
ing,  was  set  to  ten  frequency  points  for  this  study.)  The 
power  at  all  remaining  frequencies  receives  a  filter  weight 
of  one  (1.0) . 

The  processing  parameters  which  must  be  specified  for 
this  filter  are  the  low  and  high  bandpass  cutoff  frequencies, 
BPL  and  BPH,  respectively.  In  general,  these  values  were  set 
at  0.023  Hz  for  BPL,  and  at  0.059  Hz  for  BPH;  these  values 
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FILTER  WEIGHTS  OF  THE  BANDPASS  FILTER 


roughly  set  the  limits  of  the  ’signal  window’  of  the  SRO 
noise  spectra  and  exclude  the  microseisinic  noise  peak  at 
approximately  0.06-0.07  Hz  (Strauss  and  Weltman,  1977).  An 
exception  to  the  use  of  these  values  for  the  bandpass  cutoff 
frequencies  will  be  noted  in  the  disciission  of  chirp  matched 
filters . 


2.  Wiener  Filter 


The  Wiener  filter  used  in  this  long-period  signal  ex¬ 
traction  study  is  discussed  by  Lane  (1976).  Briefly,  in  the 
frequency  domain,  the  Wiener  filter  takes  the  form 


W(uj)  = 


4'ss(w) 


•tiggCu)  +  <l>„„(aj)  + 


nn 


TuT 


^ns^ 


where  W(a))  are  the  filter  weights,  (}>gg(u))  and  are  the 

autopower  spectra  of  the  signal  and  noise,  respectively,  and 
<|)g^(u))  and  are  the  crosspower  spectra  of  the  signal 

and  noise. 


It  is  generally  assumed  that  the  signal  and  noise  are 
uncorrelated  so  that  the  terms  and  (<*>)  may  be  ig¬ 

nored.  However,  there  may  be  some  correlation,  albeit  small, 
between  signal  and  noise  due  to  the  random  nature  of  the 
noise.  Attempts  by  Lane  (1976)  to  calculate  the  signal-noise 
crosspower  correlation  terms,  however,  were  not  successful 
since  the  phase  of  the  noise  is  a  random  variable.  Although 
various  models  of  the  phase  of  the  noise  were  examined  for 
use  with  the  Wiener  filter,  none  was  as  effective  as  omitting 
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the  correlation  terms  and  entirely.  The 

Wiener  filter  weights  are,  therefore,  computed  as  follows: 


W(ai) 


Two  data  items  are  necessary  for  the  successful  opera¬ 
tion  of  this  filter.  First,  a  high  SNR  event  from  the  source 
region  lander  study  is  needed  to  compute  the  (j)  (o))  term. 
Second,  a  noise  gate  immediately  preceding  the  test  signal 
gate  is  required  in  order  to  compute  the  term.  Lack 

of  high  SNR  events  or  degraded  noise  data  caused  by  instru¬ 
ment  malfunctions  will  preclude  use  of  this  filter. 


3.  Linear  Chirp  Matched  Filter 


The  signal  extraction  technique  known  as  the  linear 
chirp  matched  filter  is,  in  essence,  a  cross-correlation 
operation  of  the  form 

p(t)  =  g(T)  f(t+T)  dr 

where  f(T)  is  the  test  time  series  and  g(T)  is  a  synthetic 
dispersed  waveform. 

It  is  computationally  easier  to  apply  the  chirp  matched 
filter  in  the  frequency  domain  since  this  requires  only  a 
series  of  multiplications  rather  than  a  numerical  integra¬ 
tion.  To  see  this,  first  Fourier  transform  the  above  equa¬ 
tion: 
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PCi*>)  =  f(t+T)dTj  dt  . 

Changing  the  order  of  integration: 

P(<^)  =  .00-^”  [-00-^“  f(t+T)dt]  dT  . 

From  the  time-shifting  theorem: 

.00  -iuit  •>  j..  c/-  ^  3 w'T 

e  •’  f(t+T)  dt  =  F(u))  e-' 


and  so: 


PC^^)  =  gf"^)  P(“)  P(‘^).oo-^'*  g^'^) 


Now : 


g('^)  e'^*^"^  dt 
G(-oj)  =  g(T)  0^“"^  dx 


by  definition  of  the 
Fourier  transform 


There  fo  re 


P(a))  =  F(a))  G(-a))  . 


Since 


G(-w)  =  G*(a)),  where  G*(w)  is  the  complex  conjugate  o 
G(u)), 


P(a))  =  F(a))  G*(u)). 
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where 

is  the  lowest  frequency  of  the  chirp  matched  filter 

is  the  highest  frequency  of  the  chirp  matched  filter 

is  the  frequency  at  which  zero  phase  shift  occurs 

L  is  the  length  (in  points)  of  the  corresponding  time- 
domain  chirp  matched  filter 

At  is  the  sampling  interval  of  the  data,  in  seconds 
per  sample. 


The  frequency  domain  response  function  corresponds  to  a 
dispersive  time- domain  waveform  with  a  linear  group  delay  and 
an  essentially  flat  ampltiude  at  all  frequencies  in  the  band 
(Harley,  19  71). 

The  processing  parameters  required  for  the  operation  of 
the  chirp  matched  filter  are  the  lowest  and  highest  frequen¬ 
cies,  and  the  length  of  the  chirp  waveform. 
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4. 


Polarization  Filter 


The  polarization  filter  computes  and  applies  filter 
weights  in  the  frequency  domain  to  three -component  data  which 
have  been  rotated  to  a  vertical,  transverse,  radial  configura¬ 
tion  (see  Strauss,  1979).  The  assignment  of  filter  weights 
is  based  on  the  measured  difference  in  phase  between  the  radial 
and  vertical  components  as  compared  to  the  expected  difference 
in  phase  for  a  given  mode  of  propagation.  For  Rayleigh  waves, 
this  expected  phase  difference  is  90°.  The  filter  weights 
are  assigned  on  a  frequency-by- frequency  basis  using  overlap¬ 
ping  data  segments;  thus,  the  polarization  filter  adapts  to 
time -varying  signal  and  noise  conditions. 

A  separate  algorithm  is  invoked  in  parallel  with  the 
polarization  filter  to  extract  Love  waves  (see  Strauss,  1979). 
Briefly,  this  algorithm  determines,  on  a  frequency-by - 
frequency  basis,  the  difference  between  the  expected  arrival 
azimuth  of  the  Love  wave  and  the  actual  arrival  azimuth.  All 
energy  with  arrival  azimuth  differences  outside  some  range 
about  the  expected  Love  wave  arrival  azimuth  is  rejected  as 
noise.  Like  the  Rayleigh  wave  polarization  filter,  the  Love 
wave  polarization  filter  is  time-adaptive. 

5.  Gaussian  Narrowband  Filter 

For  the  Gaussian  filter,  the  filter  weights  form  a 
Gaussian  distribution  with  center  frequency  ,  and  with  half¬ 
amplitude  bandwidth  Aw  (see  Figure  II- 3).  The  filter  weights 
H(w)  are  defined  by: 
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FIGURE  II-3 

FILTER  WEIGHTS  OF  GAUSSIA.N  NARROWBAND  FILTER 
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HCco)  =  e  e 


+  e 


where : 


"  ■  2' 
-2a 


(W+O)^)' 

-2a 


0)^  =  center  frequency  of  the  filter 

3  =  gain  parameter;  3  =  1  to  provide  unit  gain  at  ui^ 

2 

a  =  bandwidth  parameter  related  to  the  3  dB  down 
points  by: 


+  Aw)  =0.5  H(w^)  =  0.53  ; 


thus  , 


0.53  =  3 


+  Aw/2)  -  w^)^ 


-2o 


T 


(w  +  Aw/2)  -  w  )' 
jln(0.5)  =  —5-  0 


-2a 


and 


(MI 


TOcoTT) 


Thus,  the  Gaussian  narrowband  filter  requires  two  pro¬ 
cessing  parameters:  the  center  frequency  w^,  and  the  half¬ 
amplitude  bandwidth  Aw. 


C.  THE  INTERACTIVE  SIGNAL  EXTRACTION  PROGRAM 

Implementation  of  the  interactive  long-period  signal  ex¬ 
traction  program  was  based  on  the  Extended  Interactive  Seis¬ 
mic  Processing  System  (ISPSE)  originally  formulated  by 
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Ringdal,  Shaub,  and  Black  (1975),  and  improved  by  Shaub  and 
Black  (1977).  For  this  study,  ISPSE  provides  the  analyst 
vfith  two  system  functions  and  with  the  ability  to  create  new 
user  functions  as  needed.  The  first  system  function,  SELEV, 
permits  the  analyst  to  retrieve  a  desired  waveform  from  a 
disk,  to  display  the  waveform,  and  to  select  a  processing 
time  gate.  (The  data  on  the  disk  are  obtained  from  a  mag¬ 
netic  tape  containing  waveforms  which  were  previously  edited 
and  rotated  to  a  vertical,  transverse,  radial  orientation.) 

The  second  system  function,  FILTER,  permits  the  analyst  to 
apply  bandpass  filters  and  chirp  matched  filters  to  the  data; 
to  measure  RMS  noise,  peak  amplitude,  and  surface  wave  mag¬ 
nitude;  and  to  display  the  data  after  each  operation.  Thus, 
FILTER  contains  two  of  the  five  filters  shown  in  Figure  II- 1, 
and  post- filtering  display  and  magnitude  measurement  capabil¬ 
ities  . 

To  complete  the  interactive  program,  three  user  functions 
(analogous  to  subroutines)  were  created.  These  are  POLAR,  the 
polarization  filter  algorithm,  G.^USS,  the  Gaussian  narrowband 
filter,  and  WIEN,  the  Wiener  filter.  These  user  functions 
may  be  applied  at  the  analyst's  discretion  at  any  point  after 
the  processing  time  gate  has  been  selected  by  means  of  SELEV. 
To  display  the  filtered  data  and  to  make  any  measurements  de¬ 
sired,  the  analyst  returns  to  FILTER  following  execution  of 
the  user  function. 

An  example  of  the  application  of  the  interactive  signal 
extraction  program  is  shown  in  Figure  II- 4.  To  produce  this 
figure,  the  following  steps  were  taken: 


ENSCO,  INC. 


11-12 


1 

J 


I 


r  ,  r3 


•  SELEV  -  analyst  selected  event  to  process  and 

the  processing  time  gate 

•  FILTER  -  analyst  displayed  raw  data  and  measured 

SNR  values 

•  POLAR  -  analyst  applied  polarization  filter 

•  FILTER  -  analyst  displayed  polarization  filtered 

data  and  measured  SNR  values 

analyst  applied  chirp  matched  filter  to 

polarization  filtered  data 

analyst  displayed  polarization,  chirp 

matched  filtered  data  and  measured  SNR 

values . 

An  example  of  the  application  of  the  interactive  pro¬ 
gram  to  the  problem  of  filter  parameter  optimization  is  shown 
in  Figure  II-5.  The  parameter  to  be  optimized  is  the  length 
of  the  chirp  matched  filter.  (The  sample  rate  for  all  data 
is  one  sample  every  two  seconds.)  The  top  trace  of  Figure 
II- 5  shews  the  vertical  component  of  the  Rayleigh  wave  which 
was  generated  by  a  Eurasian  event  and  which  was  recorded  at 
ANTO.  The  initial  and  final  frequencies  of  the  chirp  matched 
filter  were  set  at  0.024  and  0.077  Hz,  respectively.  In  the 
example  shown,  chirp  filters  with  lengths  of  140  ,  190  ,  240, 
and  290  points  were  successively  applied  to  the  data  shown 
at  the  top  trace,  and  the  results  were  displayed  as  shown. 

The  bottom  trace  shows  the  240-point  chirp  filter.  By  means 
of  the  parameter  measurement  module  of  FILTER,  the  SNR  values 
of  the  traces  were  then  measured.  These  results  show  that  a 
chirp  filter  length  of  240  points  produces  the  best  results 
for  this  waveform  in  the  sense  that  this  filter  length  maxi¬ 
mizes  the  SNR. 
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Processing  a  suite  o£  test  events  in  the  manner  described 
above  produced  the  optimum  chirp  filter  lengths  listed  in 
Table  II- 1  for  data  recorded  at  ANTO.  The  subregion  numbers 
used  in  this  table  are  those  defined  by  Flinn  and  Engdahl 
(1965).  The  change  in  chirp  filter  length  from  subregion  to 
subregion  indicates  the  necessity  of  making  this  type  of  fil¬ 
ter  parameter  evaluation  whenever  data  from  a  new  station- 
source  region  combination  are  to  be  processed. 

In  the  past,  attempts  have  been  made  to  determine  a 
linear  relationship  between  the  length  of  the  synthetic  chirp 
filter  and  the  epicentral  distance  (A)  of  the  test  event.  For 
example,  Strauss  and  Tolstoy  (1974)  derived  the  following 
relationship : 

CHIRP  FILTER  LENGTH  =>  10.2  7  *A-  159. 

However,  the  variance  of  the  observed  data  about  the  data  pre¬ 
dicted  by  this  relationship  was  too  large  to  permit  a  credible 
determination  of  chirp  filter  length  using  epicentral  distance. 

Table  II- 1  also  presents  the  results  of  testing  for 
optimum  chirp  filter  length  using  data  previously  processed 
by  the  polarization  filter.  Note  that  in  general,  the  opti¬ 
mum  chirp  filter  lengths  determined  for  polarization- filtered 
data  are  shorter  than  the  optimum  chirp  filter  lengths  deter¬ 
mined  for  raw  data.  This  illustrates  that  care  must  be  taken 
in  setting  the  processing  parameters  for  each  stage  of  a  cas¬ 
caded  filter  sequence  since  a  parameter  which  is  optimum  for 
one  filter  sequence  is  not  necessarily  optimum  for  another. 
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FI  inn  and  Engdahl 
(1965) 
Subregion 


Optimum 

Chirp  Filter  Length 

[n  Points 

Chirp  Filter  Ap¬ 
plied  To  Raw  Data 

Chirp  Filter  Ap¬ 
plied  To  Polariza¬ 
tion  Filtered  Data 

Rayleigh 

Love 

Rayleigh 

Love 

240 

260 

220 

260 

360 

300 

600 

500 

480 

400 

480 

260 

260 

300 

240 

240 

690 

660 

6  70 

600 

345 

335 

380 

300 
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D.  THE  AUTOMATIC  MODE  SIGNAL  EXTRACTION  PROGRAM 

The  automatic  mode  signal  extraction  program  is  struc¬ 
tured  as  shown  in  Figure  II- 1.  In  its  present  state,  the 
program  requires,  as  input,  2048-point  waveform  edits  of 
three- component ,  rotated  data.  Based  on  experience  gained 
using  the  individual  filters  and  the  interactive  signal  ex¬ 
traction  program,  the  processing  parameters  of  the  various 
filter  algorithms  can  be  set  as  follows: 

•  Two  passbands  are  sufficient  to  prefilter  the  data. 
If  a  chirp  matched  filter  is  included  in  the  cas¬ 
caded  filtering  sequence,  the  passband  is  set  at 
0.020-0.080  Hz  in  order  to  avoid  eliminating  fre¬ 
quencies  used  by  the  chirp  matched  filter.  If  the 
chirp  matched  filter  is  not  included  in  the  filter¬ 
ing  sequence,  the  passband  is  set  at  0,023-0.059  Hz 

•  The  key  input  data  item  for  the  operation  of  the 
Wiener  filter  is  the  high  SNR  waveform  required  to 
compute  <|igg(a)).  The  data  are  obtained  from  a  mag¬ 
netic  tape  which  contains  signals  with  high  SNRs 
from  each  source  region  of  interest. 

•  The  initial  and  final  frequencies  of  the  chirp 
matched  filter  are  set  at  0.024  Hz  and  0.077  Hz, 
respectively.  These  frequency  values  were  deter¬ 
mined  from  a  study  of  the  dispersion  of  signals 
recorded  at  ANTO.  The  length  of  the  chirp  matched 
filter  is  determined  from  Table  II- 1. 

•  The  polarization  filter  requires  two  parameters 
to  be  preset.  These  are  the  number  of  points  per 
processing  segment  and  the  width  of  the  pass  window 
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used  in  computing  the  filter  weights.  Optimum 
values  of  these  parameters,  determined  by  Lane 
(1976)  for  long-period  surface  waves,  were  used 
in  the  automatic  mode  version  of  signal  extrac¬ 
tion  program. 

•  The  Gaussian  narrowband  filter  requires  two  param¬ 
eters  to  be  preset.  These  are  the  center  frequency 
and  the  half-amplitude  width  of  the  filter.  The 
center  frequency  is  determined  by  the  program  by 
computing  a  smoothed  power  spectrum  of  the  data  in 
the  signal  gate  and  then  finding  the  frequency 
at  which  the  spectrum  has  its  maximum  value.  This 
frequency  is  used  as  the  center  frequency  of  the 
filter.  The  half-amplitude  width  of  the  filter  is 
set  to  one-half  the  center  frequency. 

It  is  seen  from  the  above  that  the  automatic  mode  ver¬ 
sion  of  the  cascade  filter  program  requires  minimal  input 
from  the  analyst.  For  example,  the  bandpass  filter/polariza¬ 
tion  filter/Gaussian  narrowband  filter  cascaded  filter  sequence 
(which  is  used  extensively  in  Sections  IV  and  V)  requires  only 
the  designation  of  the  signal  to  be  analyzed  since  the  filter 
parameters  are  either  preset  and  remain  constant,  or  are  de¬ 
termined  from  the  data  itself. 
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SECTION  III 

ESTIMATION  OF  SIGNAL-TO-NOISE  RATIO  IMPROVEMENT  PRODUCED 
BY  THE  APPLICATION  OF  CASCADED  FILTERS 


The  approach  used  to  estimate  SNR  improvement  produced 
by  various  cascade  filter  sequences  was  to  (1)  process  a 
suite  of  test  events  with  various  cascaded  filter  sequences; 
(2)  determine  the  SNR  gain  for  each  test  event  and  filter 
used;  and  (3)  compute  the  mean  SNR  gain  obtained  using  each 
filter  sequence.  The  test  data  selected  consisted  of  re¬ 
cordings  for  ten  earthquakes  which  were  located  in  the  source 
region  shown  in  Figure  I-l,  and  recordings  from  seven  pre¬ 
sumed  explosions  (two  from  eastern  Kazakh,  two  from  the  Ural 
mountains  area,  one  from  southwestern  Russia,  one  from  Novaya 
Zeralya,  and  one  from  central  Siberia).  All  data  were  re¬ 
corded  at  the  Ankara,  Turkey,  (ANTO)  Seismic  Research  Ob¬ 
servatory  . 


Not  all  possible  combinations  of  the  five  filters  in  the 
signal  extraction  program  were  used  in  this  test  since  some 
combinations  are  'illegal'  (illegal  combinations  are  those 
for  which  application  of  one  filter  destroys  the  signal  char¬ 
acteristic  upon  which  the  operation  of  a  subsequent  filter 
depends).  For  example,  the  narrowband  filter  cannot  be  fol¬ 
lowed  by  the  chirp  matched  filter  or  by  the  Wiener  filter 
since  the  latter  two  filters  use  the  dispersion  characteris¬ 
tic  of  the  signal  which  is  eliminated  by  the  narrowband  fil¬ 
ter. 
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The  bandpass  and  Wiener  filters  were  used  to  prefilter 
the  data;  their  primary  value  was  seen  to  be  to  cleanse  the 
data  by  suppressing  the  noise  while  minimally  affecting  the 
signal.  Using  these  filters  first  insured  that  data  pro¬ 
vided  to  succeeding  filters  were  in  the  best  possible  form. 
The  Gaussian  narrowband  filter,  on  the  other  hand,  was  used 
(with  one  exception  where  it  was  used  as  a  pre-filter  to  the 
polarization  filter)  as  a  post-filter  to  reject  all  energy 
except  that  of  the  dominant  frequency  in  the  signal  gate. 


Table  III-l  presents  the  results  of  this  SNR  gain  im 
provement  test.  The  letter  codes  in  the  column  headed 
'Cascaded  Filter  Combination'  are  defined  as  follows: 


-  0.023-0.059  Hz  bandpass  filter 

-  0.020-0.080  Hz  bandpass  filter 

W  -  Wiener  filter 

C  -  Chirp  matched  filter 

P  -  Polarization  filter 

G  -  Gaussian  narrowband  filter. 


The  order  of  the  letters  in  each  code  indicates  the  order  in 
which  the  filters  were  applied.  The  SNR  gains  of  the  indi¬ 
vidual  filters  over  the  raw  data  are  also  presented  for  pur¬ 
poses  of  comparison.  Note  that  the  Love  wave  gains  were  de¬ 
termined  only  from  earthquake  test  events  since  explosion¬ 
generated  Love  waves  are  difficult  to  detect. 


It  is  clear  from  the  data  of  Table  III-l  that  the  polar¬ 
ization  filter  dominates  the  signal  extraction  process.  In 
fact,  the  gain  of  this  filter  alone  is  greater  than  that  of 
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TABLE  III-l 

SNR  GAINS  PRODUCED  BY  VARIOUS  CASCADED  FILTER  SEQUENCES 


Cascaded  Filter 
Sequence 

LR-V  SNR  Gain 

LQ-T  SNR  Gain 

Mean 

(dB) 

Standard 

Deviation 

CdB) 

Mean 

(dB) 

Standard 

Deviation 

(dB) 

2 . 3 

1.  3 

3.2 

1.8 

BZ 

0.3 

0.5 

1.  1 

1.  3 

W 

3.  7 

2.3 

4.2 

2.6 

c 

3.4 

1.  7 

2.6 

2.5 

p 

10.5 

3.  3 

4.5 

2.7 

G 

3.  8 

2.1 

4.6 

2.6 

Bp  P 

11.4 

3.0 

4.9 

3.  8 

Bp  G 

3.9 

1.9 

5.7 

1.8 

B2,  C 

3.4 

1.  7 

2.6 

2.5 

W,  C 

4.  7 

2.5 

5.0 

2.7 

W,  P 

13.2 

4.2 

7.5 

7.0 

G,  P 

13.6 

3.  7 

5.2 

4.2 

Bi,  P,  G 

13.9 

3.  4 

6.3 

4.  3 

B2,  C,  P 

12.  3 

3. 6 

3.  0 

4.6 

B2  »  C ,  G 

5.0 

1.9 

3.4 

4.3 

B2,  P.  C 

12.2 

3.  6 

2.4 

4.0 

W,  C,  P 

13.2 

5.1 

7.  7 

4.1 

W,  C,  G 

4.4 

2.8 

5.8 

3.2 

W,  P,  C 

13.  3 

4.  7 

3.9 

7.  7 

W,  P,  G 

14. 8 

4.7 

7.  7 

6.1 

B2  >  C ,  P ,  G 

14.3 

4.  3 

5.6 

4.  1 

B2  »  P  >  C ,  G 

14.2 

4.2 

3.  5 

5.5 

W,  C,  P,  G 

14.9 

5.4 

8.  0 

5.  7 

W,  P,  C,  G 

14.9 

5.7 

2.4 

6.8 
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any  of  the  cascaded  filter  sequences  which  does  not  include 
it.  This  is  not  to  say,  however,  that  those  cascaded  filter 
sequences  which  do  not  contain  the  polarization  filter  should 
be  disregarded.  It  must  be  remembered  that  the  polarization 
filter  requires  three -component ,  rotated  (vertical,  trans¬ 
verse,  radial)  data  for  successful  operation.  Thus,  should 
the  data  be  unrotated  or  should  data  for  only  one  component 
(vertical)  be  available,  it  will  not  be  possible  to  use  any 
cascade  filter  sequence  containing  the  polarization  filter. 

In  almost  all  cases  shox\rn  in  Table  III-l,  replacement  of 
the  bandpass  filter  by  the  Wiener  filter  produces  higher  SNR 
gains.  This  increase  in  gain,  however,  is  counterbalanced  by 
the  Wiener  filter's  need  for  a  high  SNR  signal  and  a  noise 
gate  preceding  the  signal  arrival  of  the  test  event  in  order 
to  compute  the  filter  weights.  This  raises  the  point  that 
some  cascaded  filter  sequences  may  be  difficult  to  apply  in 
an  automatic  mode  of  operation  since  they  depend  on  the  avail¬ 
ability  of  information  external  to  the  test  data. 

Although  cascading  a  series  of  filters  produces  higher 
mean  SNR  gains,  the  individual  filter  gains  are  not  strictly 
additive.  This  is  illustrated  by  Table  III-2,  which  is  de¬ 
rived  from  the  data  of  Table  III-l.  In  Table  III-2,  the  'ex¬ 
pected  gain'  is  the  sum  of  the  appropriate  individual  filter 
gains  listed  in  Table  III-l,  while  the  'actual  gain'  is  the 
gain  computed  from  applying  the  cascaded  filter  sequence.  The 
letter  codes  have  the  same  meaning  as  in  Table  III-l.  This 
is  only  one  of  many  examples  which  may  be  drawn  from  Table 
III-l.  The  point  is  that  adding  more  filters  to  the  cas¬ 
caded  filter  sequence  should  not  be  expected  to  increase 
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I 
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TABLE  I I I- 2 

COMPARISON  OF  ACTUAL  AND  EXPECTED  GAINS 
OF  CASCADED  FILTERS 


Cascaded  Filter 
Stage 

Expected  Gain 

Actual  Gain 

0.3 

0.3 

®2  *  ^ 

3.  7 

3.4 

B2,  C,  P 

14.2 

12.3 

^2  >  C ,  P  ,  G 

18.  0 

14.  3 
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substantially  the  SNR  gain  of  the  sequence.  In  general,  it 
would  appear  that  the  most  efficient  cascaded  filter  se¬ 
quence  is  a  three-stage  filter. 

Based  on  the  improvement  in  SNR  gain,  on  the  degree  of 
variance  in  the  gain,  and  on  the  ease  of  application  to  a 
large,  world-wide  data  base  using  the  automatic  mode  of  op¬ 
eration,  the  favored  cascaded  filter  sequence  is  the  band¬ 
pass  filter/polarization  filter/Gaussian  narrowband  filter 
sequence  (Bj^,  P,  G  of  Table  III-l).  This  sequence  is  used 
in  Sections  IV  and  V  to  illustrate  the  effects  of  cascaded 
filtering  on  detection  and  discrimination. 

In  order  to  estimate  the  stability  of  cascaded  filter 
SNR  gains  for  different  source- receiver  paths,  the  bandpass 
filter/polarization  filter/Gaussian  narrowband  filter  se¬ 
quence  was  tested  on  earthquake- generated  surface  waves  re¬ 
corded  at  each  of  the  stations  shown  in  Figure  I-l.  The  re¬ 
sults  of  this  test  are  summarized  in  Table  III- 3  in  terms  of 
the  degree  of  RMS  noise  suppression,  the  degree  of  peak  sig¬ 
nal  suppression,  and  the  resulting  gain  in  SNR.  These  re¬ 
sults  show  that  the  SNR  gain  varies  significantly  from  sta¬ 
tion  to  station.  This  suggests  that  the  cascaded  filter  se¬ 
quence  which  yields  optimum  results  when  applied  to  data  re¬ 
corded  at  one  station  may  not  be  the  best  sequence  to  apply 
to  data  recorded  at  other  stations.  One  reason  for  this  may 
be  that  the  noise  at  different  stations  exhibits  significant¬ 
ly  different  characteristics,  and  these  effects  should  be 
taken  into  account. 
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TABLE  I II -3 

EFFECTS  OF  CASCADED  BANDPASS  FILTER/POLARIZATION  FILTER/ 
GAUSSIAN  NARROWBAND  FILTER  SEQUENCE  ON 
DATA  RECORDED  AT  NETWORK  STATIONS 


tation 

Noise 

Suppress  ion 

Signal 

Suppression 

Net 

Gain 

Mean 

(dB) 

Standard 
Deviat ion 
(dB) 

Mean 

(dB) 

St  andard 
Deviation 
(dB) 

Mean 

(dB) 

Standard 

Deviation 

(dB) 

AN  TO 

20.6 

3.4 

7.  0 

2.3 

13.6 

3.4 

GUMO 

27.3 

6.9 

12.  1 

5.9 

15.2 

7.6 

GRFO 

26.4 

4.8 

4.9 

2.4 

21.5 

5.5 

SHIO 

23.0 

5.4 

5.6 

2.9 

17.4 

3.  8 

TATO 

23.0 

4.5 

8.2 

3.0 

14.  8 

3.  5 

MAJO 

26.  8 

4.  7 

9.9 

4.9 

16.9 

7.  1 

KONO 

1 

2  3.7 

5.3 

6.0 

2.0 

17.  7 

5.0 

CHTO 

23.7 

2.8 

7.  3 

3.6 

5.3 
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SECTION  IV 
DETECTION  CAPABILITY 


A.  SINGLE  STATION  TEST 

This  section  describes  the  effects  on  single  station 
and  network  detection  capabilities  produced  by  the  applica¬ 
tion  of  a  particular  sequence  of  cascaded  filters.  This  se¬ 
quence  consists  of  a  0.023-0.059  Hz  bandpass  filter,  the 
polarization  filter  (using  the  Rayleigh  and  Love  wave  models) , 
and  the  Gaussian  narrowband  filter  (sequence  P,  G  of 
Table  III-l).  This  sequence  is  hereafter  referred  to  as  'the 
cascaded  filter.'  This  filter  was  selected  on  the  basis  of 
the  SNR  gain  produced  by  its  application  and  for  its  relative 
ease  of  application  (the  only  input  needed  was  the  designa¬ 
tion  of  each  test  event). 

Two  stations  were  selected  to  illustrate  the  single¬ 
station  case;  Ankara,  Turkey  (ANTO) ;  and  Kongsberg,  Norway 
(KONO).  The  data  base  described  in  Section  I  served  as  the 
source  of  events,  and  consisted  of  143  events  recorded  at 
ANTO  and  139  events  recorded  at  KONO.  No  presumed  explosions 
were  included  in  the  data  set  since  explosions  exhibit  a  dif¬ 
ferent  surface  wave  detectability  than  do  earthquakes.  Fur¬ 
thermore,  no  events  for  which  the  data  appeared  to  be  mixed 
or  showed  low  data  quality  (i.e.,  contained  spikes  or 
'glitches')  were  included  since  such  data  tend  to  bias  ad¬ 
versely  the  detection  thresholds. 
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In  short,  the  test  described  below  is  simply  intended 
to  demonstrate  the  improvement  in  event  detectability  pro¬ 
duced  by  the  selected  cascade  filter  sequence.  All  of  the 
detection  capability  estimates  presented  may  be  considered 
to  be  'ideal'  estimates  in  that  it  is  implicitly  assumed  in 
their  computation  that  the  analyst  will  always  see  either 
seismic  noise  or  the  sought-after  signal.  That  is,  degraded 
data  have  been  eliminated  so  that  the  analyst  is  working 
with  an  'ideal'  data  set. 

The  data  were  first  examined  in  their  raw  funfiltered) 
form  for  visually  detectable  signals  to  provide  a  detection 
baseline  against  which  the  cascade  filtered  results  could  be 
compared.  Each  event  was  judged  to  be  either  'detected'  or 
'non- detected' ,  on  the  basis  of  whether  the  data  showed 
signal-like  characteristics  (i.e.,  whether  the  data  displayed 
the  presence  of  dispersion  in  the  signal  gate,  and  whether 
the  amplitudes  in  this  gate  were  above  those  of  the  preceding 
noise  gate).  The  number  of  detected  and  non-detected  events 
for  each  bodyivave  magnitude  were  compiled,  and  the  detection 
percentage  was  computed.  Finally,  a  cumulative  Gaussian  dis¬ 
tribution  was  fitted  in  a  maximum  likelihood  sense  (Ringdal, 
1974)  to  these  detection  statistics. 

After  determining  the  detection  status  for  the  data  in 
their  raw  form,  the  cascaded  filter  was  applied  to  the  data 
for  each  event.  The  data  were  plotted  after  each  stage  of 
the  cascaded  filter  process  in  order  to  permit  the  analyst 
to  determine  the  cumulative  effect  of  each  filter  on  the 
data.  Detection  capability  curves  of  the  type  described 
above  were  then  computed  for  each  stage  of  the  cascaded  fil¬ 
ter  process . 
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The  statistic  used  to  judge  the  degree  o£  improvement  in 
the  detection  capability  is  the  50-percent  detection  thresh¬ 
old,  This  is  the  bodywave  magnitude  at  which  the  prob¬ 

ability  of  detecting  the  associated  surface  waves  is  50  per¬ 
cent.  The  analysis  results  are  summarized  in  Table  IV- 1. 

The  overall  improvement  in  detection  capability  over  that  of 
the  raw  data  produced  by  the  cascaded  filter  sequence  is  ap¬ 
proximately  0 . 8  mjj  units  for  Rayleigh  waves  and  0.6-0. 7  mj^ 
units  for  Love  waves . 

One  apparent  anomaly  in  Table  IV- 1 ,  however,  is  the  ap¬ 
proximately  0.2  mj^  units  improvement  produced  by  the  applica¬ 
tion  of  the  bandpass  filter  to  KONO  data.  One  would  expect 
improvements  of  perhaps  0.05-0.10  m^^  units  at  this  stage  (as 
indeed  is  the  case  for  the  ANTO  data)  in  light  of  the  SNR 
improvement  estimated  for  the  bandpass  filter  in  Section  III. 
Reference  to  Figure  IV- 1,  however,  resolves  this  apparent 
anomaly.  This  figure,  modified  from  data  presented  by 
Weltman,  et  al .  (1979),  shows  the  vertical  component  noise 
spectra  for  KONO  and  ANTO.  As  is  seen,  the  microseismic 
peak  of  the  noise  at  KONO  is  much  higher  than  that  at  ANTO. 
Thus,  application  of  the  bandpass  filter  used  here,  which 
rejects  the  frequencies  of  the  microseismic  peak  at  KONO, 
will  produce  higher  SNR  gains  at  KONO  than  at  ANTO.  Had  an 
interactive  system  be  used,  of  course,  the  lower  limit  of 
the  bandpass  filter  could  easily  have  been  revised  to  reject 
the  microseismic  peak  at  ANTO. 
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TABLE  IV- 1 


DETECTION  CAPABILITY  IMPROVEMENT  PRODUCED  BY  THE 
APPLICATION  OF  THE  CASCADED  FILTER  -  SINGLE  STATION  TEST 


Data  Type 

50  Percent  Detection  Threshold 
(Bodywave  Magnitude) 

AN  TO 

KONO 

LR-V 

LQ-T 

LR-V 

LQ-T 

Raw 

4.77 

4.80 

4.62 

4.75 

Bandpass  Filtered 

4.72 

4.75 

4.43 

4.57 

Polarization  Filtered 

4.12 

4.34 

3.93 

4.22 

Narrowband  Filtered 

3.99 

4.20 

3.  87 

4.05 
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0.1 


0.05 


0.025  0.0167  0.0125 

Frequency 

FIGURE  IV- 1 

VERTICAL  COMPONENT  NOISE  SPECTRA  FOR  ANTO  AND  KONO 
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B.  NETWORK  TEST 


The  approach  taken,  to  estimate  the  effects  of  the  cas¬ 
caded  filter  on  network  detection  capability  uses  a  modified 
version  of  the  network  capability  estimation  technique  re¬ 
ported  by  Snell  (1976).  This  technique  is  based  on  an  analy¬ 
sis  of  the  geometric  mean  noise  amplitudes  measured  at  the 
stations  of  the  network  (see  Table  IV-2).  The  signal  ampli¬ 
tude  at  each  of  the  network  stations  is  calculated  for  an 
assumed  event  at  epicenter  j  with  bodywave  magnitude  mj^j 
(Equation  IV-1).  With  signal  and  noise  amplitudes  assumed 
to  be  lognormally  distributed,  the  probability  that  the  SNR 
exceeds  the  individual  station  detection  threshold  is  given 
by  the  normal  cumulative  probability  function  (Equations  IV-2 
and  IV-3).  The  individual  station  detection  probabilities 
are  then  combined  into  the  network  detection  probability  of 
at  least  a  station  detections,  (>a) ,  using  a  recursive 
relationship  involving  the  individual  station  detection  prob¬ 
abilities  (Equation  IV-4).  A  detailed  description  of  this 
technique  is  given  by  Snell  (1976).  This  technique  was 
modified  so  that  instead  of  computing  detection  threshold 
magnitudes  for  the  network,  the  probability  of  detecting  the 
surface  waves  of  an  event  of  given  bodywave  magnitude  is 
plotted. 

The  geometric  mean  zero-to-peak  noise  amplitudes  for 
bandpass  filtered  data  were  taken  from  the  SRO  evaluation 
report  (Weltman,  et  al.,  1979).  To  obtain  the  corresponding 
values  for  cascade  filtered  data,  these  values  were  reduced 
by  factors  determined  from  the  mean  filter  gains  of  Table 
III-2. 
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TABLE  IV- 2 

COMPUTATION  OF  NETWORK  DETECTION  CAPABILITY 


logioCAi-)  =  m.  -  *  ^ij  (IV-1) 


P.  .  =  (j) 
ij 


logloCAi-)  -  (Pjj  +  logjpCSDT)) 


(a^  +  a^) 
^  n  S'* 


(IV-2) 


(J)(X)  = 


X  .y2/7 

f  e  y  dy 


(IV- 3) 


P.  (>a)  =  E  PJk) 
^  k=a  ^ 


CIV-4) 


where : 


b  A , c. 

A  ’  A 


Zero-to-peak  signal  amplitude  at  station  i 
for  event  j 

Bodywave  magnitude  of  event  j 

Standard  table  entries 

Station-epicenter  bias  correction 

Geometric  mean  tero-to-peak  noise  amplitude 

Variance  of  logj^O 

Variance  of  log;j^0  signal 

Normal  cumulative  probability  function 

Number  of  stations  in  the  network 

Probability  that  k  stations  will  detect 
event  j 


P. (>a)  -  Probability  that  a  or  more  stations  will 
^  detect  event  j 

SDT  -  Station  detection  threshold,  i.e.,  signal- 
to-noise  ratio  required  for  detection  at 
a  station. 
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Using  the  methods  described  above,  detection  probabil¬ 
ities  for  an  m^  =  4.0  ev'ent  with  world-wide  distribution  of 
epicenters  were  computed  in  two  ways:  first,  with  the  noise 
statistics  computed  from  bandpass  filtered  data,  and  second, 
with  the  nosie  statistics  computed  from  cascade  filtered  data. 
Figure  IV-2  shows  the  0.5  probab il ity-of-detection  contours 
for  these  two  cases.  The  dashed  lines  show  the  0.5  probability- 
of-detection  contour  derived  from  bandpass  filtered  data,  while 
the  solid  line  shows  the  0.5  prob ab il ity-of- detect  ion  contour 
derived  from  cascade  filtered  data.  The  solid  circles  (•)  de¬ 
note  the  locations  of  the  network  stations.  The  shift  of  the 
0.5  probability-of-detection  contour  away  from  the  network 
illustrates  the  effect  of  cascade  filtering  on  network  event 
detectability. 
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—  50%  Probability  of  Detection  Contour  -  Bandpass  Filter  Data  j 
; 50%  Probability  of  Detection  Contour  -  Cascade  Filter  Data  j 

5 .  ■  . 1  J; 

:rjn  60w  3ow  ow  3oe  cje  90e  i:-oe  i50e  isow 


FIGURE  IV- 2 

EFFECT  OF  CASCADED  FILTERS  ON  NETWORK  DETECTION  CAPABILITY 
WHERE  AT  LEAST  HVO  STATION  DETECTIONS  ARE  REQUIRED 
FOR  NETWORK  DETECTION  OF  AN  mj^=4.0  EVENT 
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SECTION  V 

DISCRIMINATION  CAPABILITY 


A.  SIGNAL  MEASURABILITY 

The  previous  section  discussed  the  effects  on  detection 
capability  of  the  bandpass  filter/polarization  filter/Gaus¬ 
sian  narrowband  filter  cascaded  filter  sequence  in  terms  of 
single  station  and  network  detection  capabilities.  This  sec 
tion  considers  the  measurability  of  signals  detected  after 
such  processing.  The  importance  of  examining  signal  measur¬ 
ability  is  that  the  increased  detection  capability  of  a  sta¬ 
tion  or  network  is  of  little  value  unless  reliable  magnitude 
can  be  measured  from  the  newly  detected  signals. 

The  first  step  in  assessing  the  impact  of  cascaded  fil¬ 
tering  on  the  earthquake/explosion  discrimination  problem  us 
ing  the  M^-m^^  discriminant  involves  measuring  surface  wave 
magnitudes  for  those  events  which  were  detected  after  band¬ 
pass  filtering  as  well  as  after  cascaded  filtering  using 
the  bandpass  f ilter/polarization  filter/Gaussian  narrowband 
filter  sequence.  These  pairs  of  M^  values  were  then  plotted 
as  shown  in  Figure  V-1  for  the  ANTO  data,  in  Figure  V-2  for 
the  KONO  data,  and  in  Figure  V-3  for  the  network  data  (where 
the  network  is  the  same  as  described  in  Section  IV).  A 
least-mean-square-error  fit  was  subsequently  made  to  each 
data  set,  as  shown  in  these  figures. 
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Cascade  Filter 


FIGURE  V-1 

COMPARISON’  OF  Mg  MEASURED  ON  BANDPASS  FILTERED  DATA 
TO  Mg  MEASURED  ON  CASCADED  FILTERED  DATA  (ANTO) 
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COMPARISON  OF  Mg  MEASURED  ON  BANDPASS  FILTERED  DATA 
TO  Mg  MEASURED  ON  CASCADED  FILTERED  DATA  (KONO) 
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Cascade  Filter 


FIGURE  V-3 

COMPARISON  OF  Mg  MEASURED  ON  BANDPASS  FILTERED  DATA 
TO  Mg  MEASURED  ON  CASCADED  FILTERED  DATA  (NETWORK) 
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Figures  V-1  through  V-3  show  that  surface  wave  magni¬ 
tudes  measured  after  application  of  this  cascaded  filter  se¬ 
quence  are  comparable  to  the  surface  wave  magnitudes  measured 
after  bandpass  filtering.  Since  the  least-mean-square-error 
fit  to  the  data  of  these  figures  have  slopes  of  approximately 
one  and  have  intercepts  of  0.4-0. 5,  one  can  see  that  the  ef¬ 
fect  of  cascade  filtering  on  surface  wave  magnitudes  is  to 
lower  the  surface  wave  magnitudes  by  0.4-0. 5  units.  Such 
effects  would  have  to  be  taken  into  account  in  any  applica¬ 
tion  involving  surface  wave  values. 


B.  SINGLE  STATION  DISCRIMINATION  TEST 

To  perform  this  test,  surface  wave  magnitudes  were  mea¬ 
sured  for  all  ANTQ-  and  KONO- recorded  events  which  were  de¬ 
tected  after  bandpass  filtering  and  for  all  events  which  were 
detected  after  cascaded  filtering.  (All  magnitudes  were  mea¬ 
sured  on  the  vertical  component  Rayleigh  wave.)  In  addition, 
nineteen  presumed  explosions  (listed  in  Table  V-1)  were  fil¬ 
tered,  and  surface  wave  magnitudes  were  measured  for  all  de¬ 
tections.  Inclusion  of  any  given  event  in  the  presumed  ex¬ 
plosion  list  was  based  primarily  on  its  epicenter  location; 
i.e.,  any  event  with  an  epicenter  at  a  known  test  site  was 
included  in  the  list. 

The  test  consisted  of  plotting  the  M^-mj^  pairs,  of  com¬ 
puting  a  linear  fit  to  the  earthquake  data  points,  and  of  com¬ 
puting  the  variance  of  the  data  points  about  this  linear  fit. 
The  fit  in  each  case  was  made  using  an  algorithm  which  treats 
neither  variable  as  dependent  and  which  minimizes  the  distance 
between  the  points  and  the  fitted  line  in  a  least-square  sense. 
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PRESUMED  EXPLOSIONS  EXTRACTED  FROM 
NATIONAL  EARTHQUAKE  INFORMATION  SERVICE 


Event 


Date 

Origin 

Time 

08/09/78 

18.00.00 

08/10/78 

07.59.54 

08/24/78 

18.00.02 

08/29/78 

02. 37.09 

09/15/78 

02.36.59 

09/16/78 

11.16.32 

09/20/78 

05.02.57 

09/21/78 

15.00.00 

10/08/78 

00.00.01 

10/17/78 

05.00.02 

10/17/78 

14.00.00 

10/31/78 

04.16.59 

11/04/78 

05.06.00 

11/29/78 

04 . 33.06 

12/14/78 

04.43.00 

12/18/78 

08.00 .03 

12/20/78 

04.32.58 

01/17/79 

08.00.00 

02/16/79 

04.03.56 

Latitude 

ON 

1 

Longitude 

OE 

m 

ISII 

mSM 

mKm 

67.0 

113.0 

49.0 

78.0 

50.0 

78.0 

50.0 

79.0 

50.0 

79.0 

68.0 

87.0 

62.0 

113.0 

48.0 

48.0 

65.0 

64.0 

50.0 

79.0 

50.0 

79.0 

50  .0 

79.0 

50.0 

79.0 

48.0 

48.0 

50.0 

79.0 

48.0 

48.0 

50  .0 

79.0 
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The  dashed  lines  in  each  plot  in  this  section  rep¬ 

resent  plus -or-minus  two  standard  deviations  about  the  fitted 
line.  This  measure  will  serve  to  classify  the  presumed  ex¬ 
plosions  of  Table  V- 1 ,  since  plus -or-minus  two  standard  de¬ 
viations  about  the  mean  of  a  normally-distributed  population 
encloses  95  percent  of  the  members  of  the  population.  Since 
explosions  generate  lower- magnitude  surface  waves  than  do 
earthquakes,  the  explosions  should  lie  below  the  earthquake 
population.  Therefore,  this  simple  presentation  scheme 
should  mis-classify  earthquakes  as  explosions  in  only  2.5 
percent  of  the  cases. 

Figure  V-4  presents  the  plots  determined  from 

bandpass  filtered  data  recorded  at  ANTO  and  KONO.  The  earth¬ 
quake  data  points  are  represented  by  solid  circles  (•) 

and  the  presumed  explosion  points  are  represented  by 

crosses  (+).  These  figures  show  that  only  seven  of  the 
nineteen  presumed  explosions  could  be  detected  and  classified 
based  on  the  bandpass  filtered  data  fall  seven  classify  as 
explosions).  One  earthquake  of  the  KONO  data  base  also 
would  be  classified  as  an  explosion  based  on  this  criterion. 
Note  that  only  15  percent  of  the  ANTO  earthquake  data  base 
and  21  percent  of  the  KONO  data  base  appear  in  these  plots  - 
all  other  earthquakes  were  not  detected. 

Figure  V-5  presents  the  corresponding  plots  deter¬ 

mined  from  cascade  filtered  data  recorded  at  ANTO  and  KONO. 
These  plots  show  that  fourteen  of  the  presumed  explosions  as 
recorded  at  ANTO  and  eleven  of  the  presumed  explosions  as 
recorded  at  KONO  could  be  detected  and  classified.  At  both 
stations,  however,  two  of  these  presumed  explosions  classify 
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FIGURE  V-4 

,OTS  FOR  ANTO  AND  KONO;  M  MEASURED 
ON  BANDPASS  FILTERED  DATA 


as  earthquakes.  Also,  at  both  stations,  three  o£  the  earth¬ 
quakes  would  be  classified  as  explosions.  Note  that  the  per¬ 
centage  of  the  .\NTO  data  base  appearing  in  the  plot  has  in¬ 
creased  to  52  percent  while  the  percentage  of  the  KONO  data 
base  appearing  has  increased  to  60  percent. 

Therefore,  these  data  show  that  signal  extraction 

using  a  cascade  filter  sequence  extends  the  utility  of  the 
^s""'b  tiy  increasing  the  number  of  events  which 

are  detected  and  may,  as  a  consequence,  be  classified.  Put 
another  way,  using  the  cascade  filter  sequence,  a  greater  num¬ 
ber  of  signals  would  be  detected  for  discrimination  analysis. 


c.  neuvork  discrimination  test 

To  perform  this  test,  a  subset  of  the  original  data  base 
was  first  formed  in  order  to  reduce  the  processing  time  re¬ 
quired.  A  total  of  47  earthquakes  were  selected  for  inclusion 
in  this  subset  together  with  the  nineteen  presumed  explosions. 
Each  signal  as  recorded  at  each  station  of  the  network  was 
filtered  by  the  bandpass  filter/polarization  filter/Gaussian 
narrowband  filter  cascade  sequence  and  was  examined  for  sig¬ 
nal  detections  after  application  of  the  bandpass  filter  and 
after  application  of  the  cascade  filter  sequence.  Surface 
wave  magnitudes  were  then  computed  for  each  detected  signal, 
and  the  individual  station  surface  wave  magnitudes  were 
averaged  to  form  the  network  magnitude.  Finally,  a  restric¬ 
tion  was  placed  on  the  data;  for  an  event  to  be  considered 
as  detected,  it  must  be  detected  on  the  data  recorded  by  at 
least  two  stations.  This  restriction  was  imposed  to  minimize 


ENSCO,  INC. 


V-10 


the  possibility  of  including  'surface  wave  magnitudes'  which 
were  actually  measured  on  noise. 

The  results  of  this  test  are  shown  in  Figure  V-6.  Fig¬ 
ure  V-6a  shows  the  data  measured  on  the  bandpass  filter¬ 

ed  data.  Nineteen  of  the  47  earthquakes  and  ten  of  the  nine¬ 
teen  presumed  explosions  had  a  sufficient  number  of  detections 
after  bandpass  filtering  to  appear  in  this  plot.  Note  that 
all  ten  of  the  presumed  explosions  classify  as  explosions. 

Figure  V-6b  shows  the  data  measured  on  the  cas¬ 

cade  filtered  data.  In  this  case,  44  of  the  47  earthquakes 
and  all  of  the  presumed  explosions  had  sufficient  detections 
to  appear  in  the  plot.  Fifteen  of  the  nineteen  explosions 
classify  as  explosions  while  the  remaining  four  events  fall 
into  the  earthquake  population.  Note  that  one  of  the  earth¬ 
quakes  classifies  as  an  explosion. 

Table  V-2  summarizes  the  classification  results  for  the 
nineteen  events  listed  in  Table  V- 1.  This  table  shows  that 
use  of  cascaded  filters  to  extract  signals  greatly  increases 
the  number  of  events  which  can  be  detected,  and,  therefore, 
which  are  available  for  classification  analysis. 
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TABLE  V-2 
SUMMARY  OF  DISCRIMINATION  TEST 


Event  not  detected;  no  classification  possible 
X  Event  classed  as  explosion 
E  Event  classed  as  earthquake 
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SECTION  VI 
CONCLUSIONS 


A.  SUMMARY  OF  RESULTS 

The  follow  conclusions  are  drawn  from  this  study; 

•  An  effective  long-period  signal  extraction  program 
has  been  created  and  tested  on  data  recorded  by 
Eurasian  SRO  stations.  This  program  exists  in  both 
interactive  and  automatic  mode  versions. 

•  Both  versions  of  the  program  are  modularized  for 
ease  of  insertion  of  additional  signal  extraction 
techniques . 

•  The  polarization  filter  dominates  the  cascaded  fil¬ 
ter  process.  Cascaded  filter  sequences  containing 
this  filter  yielded  SNR  gains  of  up  to  14  to  15  dB. 

•  The  bandpass  filter/polarization  filter/Gaussian 
narrowband  filter  cascade  filter  sequence  produced 
an  improvement  in  surface  wave  detection  capability 
equivalent  to  approximately  0.6-0. 8  m|^  units. 

•  Cascade  filtering  enhanced  the  discrimination 
capability  by  increasing  the  number  of  detected  sig¬ 
nals  for  which  reliable  surface  wave  magnitudes  can 
be  estimated. 
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B.  FEASIBILITY  AND  OPERABILITY  OF  APPLYING  THE  LONG- PERIOD 

SIGNAL  EXTRACTION  PROGRAM  TO  A  LARGE  WORLD-WIDE  DATA 

BASE 

Based  on  the  results  of  this  study,  the  most  feasible 
approach  to  applying  the  long-period  signal  extraction  pro¬ 
gram  to  a  large  world-wide  data  base  would  be  lo  have  both 
interactive  mode  and  automatic  mode  versions  resident  on  the 
same  computer.  This  would  permit  one  to  process  the  majority 
of  the  data  using  the  automatic  mode,  with  the  interactive 
mode  reserved  for  parameter  determination  and  the  analysis 
of  time  segments  of  more  than  usual  interest.  (It  is  anti¬ 
cipated  that  attempts  to  use  the  interactive  mode  on  all  the 
data  of  a  large  data  base  would  overwhelm  the  analyst.  Auto¬ 
matic  mode  processing  of  the  majority  of  the  data  should 
avoid  this.) 

Initially,  it  would  be  necessary  to  determine  the  signal 
suppression  of  each  selected  cascade  filter  sequence  at  each 
station  providing  data  to  the  data  base.  This  is  necessary 
to  permit  computation  of  comparable  surface  wave  magnitudes 
for  events  detected  by  different  cascade  filter  sequences, 
since  the  degree  of  signal  suppression  varies  with  the  choice 
of  cascade  filter  sequence  and  the  recording  station  (Table 
III-3) . 

C.  SUGGESTED  ADDITIONAL  RESE.i^RCH 

Significant  gains  in  detecting  and  measuring  weak  sur¬ 
face  wave  signals  can  be  achieved  by  evaluating  and  optimizing 
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an  alternate  design  for  the  cascade  filter  used  in  the  signal 
editing  process.  The  cascaded  bandpass  filter/polarization 
filter/Gaussian  narrowband  filter  sequence  (B^ ,  P,  G;  see 
Table  III-l)  evaluated  in  this  study  is  considered  the  opti¬ 
mum  choice  of  a  processor  which  automatically  edits  and  mea¬ 
sures  surface  wave  signals.  It  was  designed  for  maximum- 
likelihood  detection  of  at  least  85%  of  the  signals  which  are 
expected  to  be  observed. 

Another  cascade  configuration,  the  Wiener  filter/polar¬ 
ization  filter/chirp  filter/Gaussian  narrowband  filter  se¬ 
quence  (W,  P,  C,  G] ,  provides  an  enhanced  capability  for  ex¬ 
tracting  weak,  hard-to-detect  signals  which  are  missed  by 
the  P,  G  sequence.  The  capability  of  the  W,  P,  C,  G  pro¬ 

cessor  to  detect  weak  signals  at  the  15%  level  of  detection 
results  from  the  higher  variance  associated  with  this  proces¬ 
sor.  This  is  especially  the  case  for  the  network- related 
performance  of  the  W,  P,  C,  G  processor  in  that  at  least  one 
additional  detection  out  of  the  eight  SRO  stations  used  here 
is  expected  at  the  15%  level  of  detection. 

For  enhanced  network  performance,  it  is  proposed  that 
the  W,  P,  C,  G  cascade  filter  sequence  be  operated  in  paral¬ 
lel  with  the  Bj^,  P,  G  cascade  filter  sequence.  Based  on  the 
statistics  of  Table  III-l,  it  is  expected  that  at  least  0.17 
units  of  improved  magnitude  detection  capability  can  be 
achieved  in  this  way  for  single-station  detections.  Consid¬ 
erably  more  gain  is  expected  to  be  achieved  when  parallel 
processors  are  used  with  a  network.  This  can  be  accomplished 
(with  only  a  minor  increase  of  computer  time)  by  utilizing  the 
W,  P,  C,  G  cascade  filter  sequence  as  a  backup  to  the  automated 
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,  P,  G  cascade  filter  sequence.  In  this  dual  mode  of  opera¬ 
tion  the  W,  P,  C,  G  cascade  filter  sequence  would  be  used 
only  when  signals  are  not  detected  after  application  of  the 
P,  G  cascade  filter  sequence. 
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